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Abstract 

The dd — > 3 Hen reaction is considered at the energies between 200 MeV and 520 
MeV. The Alt-Grassberger-Sandhas equations are iterated up to the lowest order 
terms over the nucleon-nucleon t-matrix. The parameterized 3 He wave function 
including five components is used. The angular dependence of the differential cross 
section and energy dependence of tensor analyzing power T20 at the zero scattering 
^D ■ angle are presented in comparison with the experimental data. 
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1 Introduction 

During several decades hadronic reactions with helium and tritium were extensively inves- 
tigated at the energies of a few hundred MeV. A number of experiments to study nucleon 
knockout from polarized 3 He were performed at TRIUMF. As a result, differential cross 
sections were measured at the energies between 220 and 590 MeV p],[2],[3]. Moreover, the 
polarization observables, such as analyzing powers and spin correlation parameter, were 
obtained at 220 MeV [2] and 290 MeV pQ. The analyzing powers and spin correlations 
were also studied at IUCF at the energy of 197 MeV [5]. 

The aim of these experiments was to study the helium internal structure. The simple 
relations between the helium wave function and differential cross section and polarization 
observables in the frame of the plane- wave- impulse-approximation (PWIA), give an op- 
portunity to extract useful information about the ground state spin structure of helium. 
In order to study the high-momentum components of the 3 H e, the elastic backward scat- 
tering of p 3 He was investigated at RCNP (Osaka). Here the differential cross section 
P^ ! and spin correlation parameter C yy were measured at proton energies of 200, 300, and 
400 MeV [6]. 

Several years ago an experiment to study dd — > 3 Hen, dd — > 3 Hp reactions was carried 
out at RIKEN [7J, [8]. The vector and tensor analyzing powers were obtained in a wide 
angular range at three deuteron kinetic energies: 140, 200, and 270 MeV. Previously the 
differential cross sections of the reactions 2 H(d, n) 3 He and 2 H(d, p) 3 H were measured in a 
wide angular range for incident deuteron momenta between 1.1 GeV/c and 2.5 GeV/c [9]. 
The dd — > 3 Hen reaction was considered in the one-nucleon-exchange (ONE) framework 
in ref. [iTRfTT]. High sensitivity of some of the polarization observables was shown to the 
spin structure of the 3 H e. However, the data obtained at RIKEN are in disagreement 
with ONE predictions. Only a small angular range, around 0° — 15° and 165° — 180°, 
is reasonably described by ONE mechanism. This result stimulated further theoretical 
investigations of this reaction. 

The four- nucleon problem is topical up to now in spite of many efforts to solve it. 
Significant progress in the studies of the dd — > 3 Hen(tp) reactions was achieved at low 
energies (Ed < 5 MeV) [12], [13], [14] . Here the reasonable description of the experi- 
mental data was obtained both for the differential cross sections and for the polarization 
observables. 



Practical integral equations for the four-body scattering were developed by Grass- 
berger and Sandhas [T5] • In this formulation the original operator relations were reduced 
to effective two-body equations in two steps by employing separable expansions both for 
the two-body and for the three-body subamplitudes. After the partial wave decomposition 
we deal with one-dimensional equations. 

This approach was applied in ref. [16J to study dd — > pt{n 3 He) reactions and p 3 He 
elastic scattering at the energies up to 51.5 MeV. Here the first-order i^-matrix approxi- 
mation was applied to solve effective two-body equations. The obtained results reasonably 
describe the shape of the differential cross sections but fail to reproduce the second max- 
imum in the differential cross section of the dd — > pt(n 3 He). Inclusion of the principal 
value part of the propagators and use of different potentials did not result in significant 
improvement p2] . Nevertheless, the carried out investigations have shown that the agree- 
ment between the theoretical predictions and data improves with increasing the energy 
when the second maximum is not so evident. 

At higher energy the four-nucleon problem was considered in ref. [18] , where the deuteron- 
deuteron elastic scattering was studied at 231.8 MeV. The approximation based on the 
lowest order terms in the Neumann series expansion of the AGS- equations, was used 
to describe the differential cross section and vector and tensor analyzing powers. The 
obtained results have demonstrated the underestimation of the differential cross section 
while the curves for the deuteron analyzing powers reproduce the behaviour of the data 
at forward angles. 

In the present paper the dd —¥ 3 Hen reaction is studied at the deuteron energies 
between 200 MeV and 520 MeV. We start our investigation from AGS equations for the 
four-body case [15] and then iterate them up to the first order terms over the nucleon- 
nucleon t-matrix. In such a way we include not only ONE mechanism into consideration 
but also the next term. It corresponds to the case when nucleons from different deuterons 
interact with each other and then form a three-nucleon bounded state and a free nucleon. 
The parameterization based on the modern phase-shift analysis data is applied to describe 
NN interaction. The partial wave decomposition is not used in this approach. It allows 
us to avoid the problem related with convergence which is important at the considered 
energies. 

The paper is organized as follows. Section 2 gives the general formalism. Here the 
expansion of the AGS equations is presented for the dd — > 3 Hen reaction. In this section 
the 3 He wave function is discussed, and the description of the nucleon-nucleon interaction 
is presented. The details of calculations of the scattering amplitude terms are also given. 
The obtained results are discussed in Sect. 3. The conclusions are contained in Sect. 4. 

2 General formalism 

Here we consider the reaction where four initial nucleons are bounded in pairs forming 
two deuterons, and three final nucleons are bounded to the helium or tritium and one 
nucleon is free. In other words, we have the reaction of the (2) + (2) — >• (3) + (1) type. 
We write the transition operator U(z) for our reaction as it was offered by Grassberger 



and Sandhas [151: 



Uf3 a (z) = (1 - 6j3 a )(z - H a ) + ^J T ik (z)G (z)Uik, a (z) + ^ V^,^ , (1) 

where a and /3 denote two-cluster partitions of the four-particles. Here these labels are 
referred to the initial and final states, respectively: 

l)a = (ij)(kl); V a = V^^i) = V%j + V M ; $(y)(fcz) = \kj, hi > \4>ij > \i>ki > (2) 
2)/3 = (ijk); Vp = % fc) = Vij + V jk + V ik ; $ fo - fc) = \k ijk , k t > \tp ijk > . 

In accordance with the AGS-formalism the channel Hamiltonian is defined as a sum of 
the free particles Hamiltonian Hq and the interaction potential: 

H a (/3) — H + V a (i3). (3) 

The eigenf unctions of the channel Hamiltonian |$ a > characterize possible initial 
and final configurations. These functions are products of plane waves and internal wave 
functions \ip a >. 

The operator Tij(z) in Eq.(TfJ) is a two-body transition operator which satisfies the 
Lippmann- Schwinger equation: 

T ij (z) = V ij + V ii Go(z)T ij , (4) 

where Go is the resolvent of the four-nucleon kinetic energy operator Gq(z) = (z — Hq)^ 1 . 

The operator Uik, a in Eq. (CD) corresponds to the case when the initial state a is deter- 
mined as in Eq.(j2]) and the final state is a combination of two bounded nucleons (ik) and 
two free nucleons. This transition operator can be also defined from Eq. ([1]) if we put the 
final state (3 = (ik). The notation ik ^ P means that pair (ik) is not either equal to one 
cluster of /3 or contained in it. 

We deal with four identical nucleons and two identical deuterons in the initial state. 
It means that symmetrized wave functions both for the initial and final states, should be 
built. Following ref. [19] we have constructed a wave function for the initial state where 
four nucleons form two bounded states: 

1^(12)^(23) > s = IJjL£qM12) >s 1-0(23) > s . (5) 

Q 

By Q we denote all possible permutations of two nucleons. For four particles we have 4! 
permutations of this kind that is reflected in the first factor. The second coefficient is 
from normalization of the symmetrized wave function. The wave functions of deuterons 
i/j(ij) s are also antisymmetrized: 

\iP(ij) > s = A=MiJ) > -Mi) >]• (6) 

Three nucleons in the final state are bounded and one nucleon is free. The corresponding 
symmetrized wave function is as follows: 

1-0(123)4 > s = ^ J2 Q\^( 123 ) >s K > • (7) 



Here three-nucleon state (ijk) s is also presented by the antisymmetrized wave function 

|-0(123) > s = -=[\ip(123) > -1-0(213) > H-|-0(231) > —[-0(321) > +\ip(312) > — 10(132) >]. (8~ 
v6 

After straightforward calculations the reaction amplitude can be written as: 
< n 3 He\U\dd >= ±=[< 4,V(123) S |C/|V(12) S V(34) S > - < 1, 0(234) S |C/|0(12) S 0(34) S >]. (9) 

The same way it is necessary to find two matrix elements of the transition operator U. We 
start to consider the first of them. This term corresponds to the case of /3 = (ijk) = (123), 
a = (12) (34). From Eq.flD we get 

^(123),(12)(34) (z) = {Z - H ) ~ Vn - V 34 + T U (z)G (z)U( U ^i2)(34:) (z) + 

T 24 {z)G (z)U(24),(12)(34){z) + T 34 (z)G (z) C/(34),(12)(34) • (10) 

This relation contains transition operators for another reaction type. In the final state 
two particles are bounded and the other two are free, while the initial state is the same 
as before. In order to derive expressions for these operators, it is convenient to rewrite 
Eq.(TjQ) in the following form: 

U/3a{z) = (1 - 8p a )(z - Hp) + ^ Up, mn {z)G (z)T mn (z) + ^ V P$l3,mn- (H) 

Then by putting (3 = (ij) we obtain 

^(U),(i2)(34) {z) = {z- H ) + U {14)i{13) (z)G (z)T 13 (z) + U {14)t{14) (z)G (z)T 14 (z) + 

U(u),(23){z)G {z)T 2 z{z) + U(u),(24){z)Go{z)T 24 {z) 
U(24),(i2)(3<±)(z) = {z- Ho) + U(24),(i3)(z)Go(z)T 13 (z) + U( 2 4),(u){z)Go{z)T u (z) + (12) 

U(24),(23){z)G (z)T 2Z {z) + U {24 ),(24){z)G (z)T 24 (z) 
U(34),{12)(34){z) = {Z- Ho) - V 34 + U(34),(13){z)G {z)T 13 (z) + J7(34),(14) {z)G {z)T 14 (z) + 

U(u),(23)(z)Go(z)T 23 (z) + U {34)t(2 4)(z)Go(z)T2 4 (z). 

Iterating these equations only up to the first order term over T-matrix, we get the following 
sequence for the £/(i23),(i2)(34)-operator: 

^(i23),(i2)(34) « {z - Hu) + T l4 {z) + T 24 (z) . (13) 

Likewise we derive the expression for the other transition operator in Eq.( fT9l) : 

*/(234),(12)(34) « {Z ~ H 34 ) + T l3 { Z ) + T u (z). (14) 

Since the initial and final states are antisymmetrized, the contributions of the T 24 and 
T 14 matrix elements are equal to each other. In order to show it, we use the properties of 
the permutation operator: P12-P12 = 1, P12P24P12 = T 14 . 

< 4, V(123) S |T 24 |V(12) S V(34) S >=< 4, ^(123) S |P 12 P 12 T 24 P 12 P 12 |^(12) S ^(34) S >= 
= < 4, ^(213) S |T 12 |^(21) S ^(34) S >=< 4, 0(123) S |T 12 |0(12) S 0(34) S > (15) 



It also concerns T 13 and T 14 matrix elements in the exchange contribution: 

<l,V(234),|Ti 3 |^(12)W>(34),> = < l,^(234) s |Ti 4 |^(12)^(34) s > . (16) 

On the other hand, using the permutation operator P 14 we can get the following useful 
relation: 

< l,^(234),|T 14 |V(12)W(34) a > = - < 4,^(123) s |P 14 T 14 |t/;(12)^(34) s > . (17) 

It gives us an opportunity to join all terms with NN T-matrix into one. In such a way 
Eq. (fl9l can be reduced to the following: 

< n 3 He\U\dd >= -=[< 4, ip{123) a \z - H 12 \^(12) s ip(34) a > - (18) 

V6 

< 1,^(234)^-P 3 4^(12)^(34) S > +2 <4,^(123) s |T 1 s r|^(12) s ^(34) s >], 

where antisymmetrized NN T-matrix is defined as Ty[ m = (1 — Pi 4 )Ti 4 . 

In order to simplify the statement below, we divide the latter expression via three 
terms corresponding to these contributions into the reaction amplitude: 

<n 3 He\U\dd> = ^6(2E d -E h -E n ){J° r NE + J°» E + 2J ss }. (19) 

All the calculations have been performed in the center-of-mass. The following definitions 
are introduced for the momenta and energies of the deuterons, helium, and neutron: 

P i= -p 2 = p d , P h = -p n (20) 



E d =^M* + Pl E h = ^M£ + P£, E n = ^mf,+pl 

Two first terms in Eqs. (TT8l) . (TT9|) correspond to the one-nucleon-exchange (ONE) mecha- 
nism of the reaction. We call the first of them as "direct" and the second one as "ex- 
change" . Here one of the deuterons breaks in a neutron and proton. One of the nucleons 
becomes free, while the other interacts with the remained deuteron forming helium or 
tritium. Schematically it can be presented by diagrams in Figs. la and lb. The latter 
term corresponds to single scattering (SS) when two nucleons from different deuterons 
interact in the final state (Fig, lc). 

2.1 One-nucleon-exchange 

We start from consideration of ONE contributions. Taking the quantum numbers and 
momenta of all particles into account, we get the following expression for ONE terms: 



<n 3 He\U\dd> = —<p n m n T n \ < * 123 (P h M h r h )\(2E d - J ff 12 )|$ 12 (P 1 ,M 1 ) > |$ 34 (P 2 ,M 2 ) > - 

<p n m n r n \ < V 2U (P h M h T h )\(2E d - H 34 )\^ 12 (P 1 ,M 1 ) > |$ 34 (P 2 ,M 2 ) > . (21) 

Here we introduce notation \l/^ fc (P/j, Mh, r^) for the 3 He wave function, where 3 He 
is formed by i,j,k nucleons and has momentum P^, spin projection Mh and isospin 





(a) (b) 




(c) 

Figure 1: The diagrams taken into consideration: one-nucleon-exchange (a),(b), and single 
scattering (c) graphs. 

. . — * 

projection Th- Note in case r = — 1/2 we deal with the reaction of dd — > tp. The $* J (P, M) 
denotes the wave function of the deuteron with momentum P and spin projection M . 
Inserting the unit operator into Eq.( l2T|) : 



1 = J dpmdpW^r u p 2m2 r 2 ,p 3m3 r 3 X pW^-Ws^l , (22) 

for ONE contribution we get 

< n 3 He\U\dd > = — p / dpidp 2 dp 3 

(E d - E n - E 3 ) < y 123 (P h M h T h )\p 1 m 1 T 1 ;p 2 m 2 T2;p 3 m 3 T 3 > 
<p 1 m 1 T 1 ;p 2 m 2 T 2 \<5> 12 (P d ,M 1 ) >< p 3 m 3 T 3 ;p n m n T n \<5> u {-P d ,M 2 ) > - 

- — = / dp 2 dp 3 dp 4 (23) 

(E d - E n - E 2 ) < ^ 2M (PhM h T h )\p 2 m 2 T 2 -,p 3 m 3 T 3 -,p 4 m 4 T 4 > 

< Pnm n T n ;p 2 m 2 T 2 \<5> l2 (P d , M x ) >< p 3 m 3 T 3 ;p 4 m 4 T 4 \<5> u {-P d , M 2 ) > . 

Henceforth, we imply summations over all dummy discrete indices. 

In our calculation we use the parameterized wave function of a three-nucleon system 
offered in ref.[20j. This wave function was derived by fitting the full Faddeev wave function 
obtained with the CD Bonn [2T] and Paris [22] NN-potentials. The wave function is fully 
antisymmetrized and defined in terms of the nucleon pair and spectator momenta. If 



Table 1: Quantum numbers of the partial waves included into the definition of the three- 
body wave function. T, S, C, J refer to the isospin, spin, orbital momentum and total 
angular momentum of the NN-subsystem. I is the relative orbital momentum of the 
spectator and K is the channel spin [20] . 
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we choose particles (12) as a pair and particle 3 as a spectator, the three-nucleon wave 
function is presented in the following form: 

< p qu\^ > = < p qV\i/j[(12)3] > + < jwIpWi^s >< P23gi^2 3 |^[(23)l] > +(24) 

+ < P qv\P31 92^31 >< P31?2^3l|V ; [(31)2] > . 

Here the following notations have been introduced for pair relative momentum p Q and 
spectator momentum q 

Pi — P2 P -> 

Q = Pa - "g , P = Pi+P2+P3 = Ph 

qi=Pi-j (25) 

_ P 

Q2 = Pi - "g" • 

The radial part of the three-nucleon wave function (I24"j) is presented as a sum of the 
two terms each of them has a separable form: 

< poH^ >= «i (PoK(g) + *£(poK(?), (26) 

where ?; A (po), "^a are defined as: 

^Saw- *£?w' x -^- (27) 

Index i/ denotes number of one of the three-nucleon channels (Table 1). The five chan- 
nels are included into the definition of the wave function: 1 s S, 3 SiS, 3 SiD, 3 d 1 S, 3 diD. 
Parameters a v n A , b v n A , m v n A and M% A can be found in [20J. 

The wave function of the deuteron, which contains (ij) nucleons, is denoted in Eq. (j2T|) 
as <&(P d ,M d ). Here P d and M d are momentum and spin projection of the deuteron, 
respectively. In the rest frame the non-relativistic wave function of the deuteron depends 
only on one variable po which is the relative momentum of the proton- neutron pair: 

< P d /2 + p , m p , P d /2 - p m n \$(P d M d ) >= (28) 

J2 u l(Po) < ^rn p -m n \\M s >< LM L lM s \lM d > Y^ L (p ) , 

L=0,2 



where u (p ) and u 2 (po) describe S and D components of the deuteron wave function [2T] . 
[23], Po is the unit vector in p direction, and m p , m n are the proton and neutron spin 
projections, respectively. 

Using transformations of vectors pi, p 2 , f>3 to Jacobi variables ( 125]) and taking into ac- 
count momentum conservations in the deuterons and helium vertices, we get the following 
expression for the first term in Eq. (l2TT) : 



JZ NE = < 1 -M h 1 -m n \T^\lM 1 lM 2 >- 



■1) 



1/2-Tft 



V2 



K, j dpodpoplripo, q)Yt\q)Y^\p Q ) < CM1M s \JMj > (29) 

< JMj-m 3 \KM K >< l\xKM K \~M h >< -m 3 -m n \lM' s >< LM L \M S \IM 1 > 

< L'M' L 1M' S \1M 2 > u L (p )Y^(p )u L ,(\P h - P d /2\)Y^(P h ~^~P d /2) 
with kinematical factor /C defined as 



K = E d - E n - \]m N + (P h - P d y. (30) 

Definitions Eq. (j2"3]) and (|2"8"]I have been also used to obtain this equation. Superscribe 
index v of the helium wave function marks one of the five channels considered in J2U] arid 
defined by quantum numbers of the nucleon pair (£, J), the relative orbital momentum 
of the spectator I and the channel spin K [21] (Table 1). We also preserve here the 
dependence on isotopic number r^ that allows us to consider both dd — > 3 Hen and dd — > tp 
reactions. As it follows from Eq. (l2"5"]) . spectator momentum <f is defined only by helium 
and deuteron momenta, q — \Ph ~ P d . Since only two spherical functions in Eq. (!29|) are 
dependent of integration angles, we can simply integrate this expression over the angular 
dependence of p : 

J° r NE = L ^ K j d PoPlr(P0,q)u L {p )u L/ (\P h - P d /2\) < lM^mzlKMK&l) 

< lfiKM K \^M h >< L'M' L 1M' S \1M 2 >< l -m 3 l -m n \lM' s > Y^(P h ^~P d /2)Yr(q) ■ 



After substitution of the partial wave decomposition of the helium wave function [20] . 
we get the final expression for the direct term of the ONE-contribution: 



rONE __ {-1) 1/2 ~ T V f^^M\Ph-PdM) . 1_ 1. 



jy,r = v ^ k. j dp pi r^^ a '^ < -m 3 - mn \iM 2 > + ( 32 ) 



u 2 (\P h - P d /2\)Y 2 ML (P h - P d /2) < 2M' L 1M' S \1M 2 >< -m 3 -m n \lM' s >} 

1 11 

{—?=[uo{pom{Po, q) + u 2 (Po)^Pi(,Po,q)} < lM l7 -m 3 \-M h > + 

V47T ^ £ 



[u (p )^(po,q)+u 2 ip )^(po,q)}Y^ *(q) < lM 1 -m 3 \-M K >< 2^-M K \-M h >} . 

This expression contains only four components of the helium wave function, since channel 
v = 1 corresponds to the isotriplet state of the pair which is forbidden for the ONE- 
mechanism. 

8 



In order to get the exchange term of the ONE-amplitude J°Jch -, it is necessary to 
replace Pd — > —Pd and M 1 < — > M 2 in the previous expression. 

2.2 Single scattering 

The single-scattering term in Eg. (1181) can be rewritten in a more evident form: 

J ss = dp 1 dp 2 dp 3 dp 4 d]? 1 <y l23 {P h M h T h )\]? 1 m[T[]p 2 m 2 T2]P3m 3 T 3 > 

< p , 1 m[r[;p n m n T n \T(2E d - E 2 - E^l^m^; p 4 m 4 r 4 > (33) 

<PimiTi;p 2 m 2 T2|<l> 12 (Pd,Mi) ><P3m 3 T 3 ;p 4 m 4 T 4 \^ M (-P d ,M 2 ) > . 

We have here five integration vectors but three of them can be removed due to the momen- 
tum conservation. We introduce vectors ko and k' which correspond to the neutron-proton 
relative momenta in the deuterons: 

£o = g (Pi ~ Pa)> ^o = j (Pa ~ P^- ( 34 ) 

As it is mentioned above, we have used the three-nucleon wave function in a separable form 
which depends on two variables: Pq, a relative momentum of a pair, and spectator momen- 
tum q. In our calculations it is convenient to choose the nucleon pair (23) as a cluster and 
nucleon 1 as a spectator. It is possible since our wave function is symmetrized:^ 123 = \1/ 231 . 
Then the arguments of the helium wave function are expressed via momenta ko and k' : 

1 -. 2 -> 

Vo = ^(Pd-k -k' ) q = -P h + k -k' . (35) 

Using the definitions of 3 He and deuteron wave functions, Eqs.(l24l).( |28l) . we can write 
the following expression for the SS-term: 

jss =< l Mh ^m n \T ss \lM 1 lM 2 >= 

( ~ 1) 2 T1 ^ J dk dk' o r(po,q)Yf(qyY^(p o y < lfiKM K \^M h > 

1 11 

< CMcSM s \ JMj >< JMj-m^KMx >< -m 2 -m 3 \SM s > 

< TM T l -T[\^r h ><\~ n^T 3 \TM T >< ±r^ - r h \T'M' T >< X -r^- - r 3 \T ' M' T > 
<P h + k - k' , m[; -P h , m n \T(2E d - E 2 - E 3 )\P d /2 + k , m i; -P d /2 - k' , m 4 > 

< -m x -m 2 \lM >< LM L 1M\IM 1 > u L (k )Y^ L (k ) (36) 

< -m 3 -m 4 \\M' >< L'M' L IM'\IM 2 > u L/ (k' )Y^ L (k' ) . 

The nucleon-nucleon scattering is described by the T- matrix element. We use the 
parameterization of this matrix offered by Love and Franey [25]. This is the on-shell NN 
T-matrix defined in the center-of-mass: 

< KVl/4l*c.m.|«>lA<2 > = < kVi/41^ 4 + 5(<T!iV*)(cr 2 iV*) + (37) 

C(<n + <r 2 ) ■ N* + DitnDfaq*) + F(mQ*)(«r a Q*)|K>i^2 > • 



The orthonormal basis {g*, Q*, N*} is a combination of the nucleon relative momenta in 
the initial k* and final «/* states: 

K,* — K*' - K* + «*' - K* X K*' , , 

<f = i r> Q = i — I — n, N = i r- ( 38 ) 

The amplitudes A,B,C,D,F are the functions of the center-of-mass energy and scattering 
angle. The radial parts of these amplitudes are taken as a sum of Yukawa terms. A new 
fit of the model parameters [26J was done in accordance with the phase- shift- analysis data 
SP07 [27]. 

Since the matrix elements are expressed via the effective iViV-interaction operators 
sandwiched between the initial and final plane-wave states, this construction can be ex- 
tended to the off-shell case allowing the initial and final states to get the current values 
of k and k'. Obviously, this extrapolation does not change the general spin structure. 

In order to relate c.m.s. and the frame of our calculations, first of all, we apply Lorentz 
transformations to kinematical variables. Let us consider momenta and energies of the 
colliding nucleons: 

Vi = k* + u(^1 + E*Y E t = jE* + (uK*) (39) 

p 4 = -k* + u C " x +E*), E A = 1 E* ~{uk*), 

where E* is the energy of one of the nucleons in c.m.s. By u = (7, u) we denote the 
4-velocity of the reference frame relatively c.m.s: 

_ Pi+Pa E x + E± 

u = -=-, 7 = -=—. (40) 

V s V s 

Mandelstam variable s is defined as usual: 



V T S = 2E* = y/(E 1 + £ 4 ) 2 - (ho - Kf (41) 

Then two-nucleon state in the reference frame can be related with that in the c.m.s. due 
to rotations in the spin space of these nucleons: 

\p 1 m 1 > |p 4 m 4 >= N\n*m' x > I - K*m' A > £ m ' imi («, j?i).D m ^ m4 (tt,p4), (42) 

where the Wigner rotation operator is 

_ u ( a[u X p\ \ 

D(u,p) = expi^au} = cos- (l - , (1 - -j- - - - ^ . (43) 

The rotation is performed around the axis [u x p\ on the angle u determined by 

tg-z = 7T- /' r 44 

2 (1 +7)111 + E+ y/S 

After transformations (f39l the two-nucleon relative momentum in c.m.s. n* is written as 
follows: 

K * = Pd W + V~s) + k' Q (E 1 + yj/2) 

2 E 1 + E A + y /s ' l ' 
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Figure 2: The differential cross section at the deuteron momentum of 1.109 GeV/c as a 
function of t. The data are taken from [9]. 

Likewise we can obtain an expression for the relative momentum of the scattered nucleon 
pair: 



K*> = P h + 

Here we use the following definitions: 



(k -k' )(E n + ^s'/2) 



(46) 



Ei = \Jm> + (fc + P d /2)\ E, = ^m2 + (k' + P d /2)*, 

E[ ---- Jm* + (PI + k - k' )\ V^ = 2E'* = J(E[ + E n f - (k - ktf 



(47) 



It is also necessary to take into account normalization factor M and a kinematical 
coefficient due to the transformation of the off-energy-shell T-matrix [28] . The expressions 
for these factors were given in detail in ref.|29j. 

It should be noted, that the current parameterization describes the NN-interaction in 
the wide energy range between 50 MeV and 1100 MeV [26J. However, at low energies 
(< 100 MeV) the quality of the parameterization can not be assessed due to the lack of 
the experimental data. Therefore, we do not consider the single scattering contribution 
into dd —¥ 3 Hen reaction amplitude at the deuteron energy below 200 MeV, where the 
Faddeev calculation technique is more preferable. 
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Figure 3: The differential cross section at the deuteron momentum of 1.387 GeV/c as a 
function of t. The data are taken from [9]. 

3 Results and discussions 

The formalism presented above was applied to describe the experimental data obtained 
for dd — > 3 Hen and dd — > tp reactions at the deuteron kinetic energies of a few hundred 
MeV. The calculations have been performed with CD-Bonn deuteron and helium wave 
functions. The differential cross section can be written as a function of Mandelstam 
variables s and t = (Pd — Phb) 2 '- 



da 4vr s 2 - {M 2 He - m 



dt 9s 



2\2 



s - AM} 



■\J(s,t)\" 



(48) 



We consider three energies, 300 MeV, 457 MeV and 520 MeV, which correspond to 
the laboratory momenta P^ = 1.109, 1.387, and 1.493 GeV/c, respectively. In this 
energy range the presented formalism is more successful. Moreover, we have a set of the 
experimental data on the differential cross sections in a wide angular range obtained at 
these energies in Saclay [9]. 

In Figs. 2-4 the results of the calculations of the differential cross sections are presented 
in comparison with the data. In order to demonstrate the contribution of the single scat- 
tering term, we have considered two cases. One of them corresponds to the calculations 
including only ONE terms. The results of these calculations are given with the dashed 
curves. The other case corresponds to the calculations taking into account both ONE and 
single scattering contributions. These results are presented with the solid curves. 
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Figure 4: The differential cross section at the deuteron momentum of 1.493 GeV/c as a 
function of t. The data are taken from [9]. 

As expected, the contribution of the rescattering term is not large at small scattering 
angles (t ~ 0.7 — 0.8(GeV/c) 2 ). It is in agreement with the results obtained in ref. [8]. 
However, the difference between these two curves increases with the angle and reaches the 
maximal value at 90°. Taking the single scattering diagram into consideration significantly 
improves the agreement between the experimental data and theoretical predictions. We 
have a good description of the data for Pi a b = 1.109 GeV/c (Fig. 2). Nevertheless, the 
underestimation of the differential cross sections is observed at the deuteron energies 
above 300 MeV (Figs. 3, 4). Perhaps, this discrepancy can be reduced, if the A-excitation 
in the intermediate state is taken into account. This possibility is discussed in ref. [9], 
where the A-isobar is taken into consideration in the simplest phenomenological model. 

The formalism presented here gives us an opportunity to calculate not only the differ- 
ential cross sections but also polarization observables. In this paper we have considered 
the energy dependence of tensor analyzing power T20 at the scattering angle equal to zero 
(Fig. 4). The experimental data were obtained at RIKEN [7j. As it is mentioned above, the 
contribution of the single scattering term is not large at small angles. Nevertheless, one 
can observe some improvement of the agreement between the data and theory predictions. 
Unfortunately, we do not have enough experimental data to confirm this tendency. 

4 Conclusions. 

The model to describe the dd — > 3 He n(tp) reaction at the energies of a few hundred 
MeV has been presented in this paper. We start from the AGS-equations for N-body 
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Figure 5: The energy dependence of tensor analyzing power T 20 at the zero scattering 
angle. The data are taken from [7] 

system. Iterating these equations over the NN t-matrix we obtain the expression for the 
reaction amplitude. In the presented calculations only two lowest terms of this expansion 
are included into this consideration. Here we do not solve any equations to define wave 
functions of the bounded states or to find the nucleon-nucleon t-matrix. Instead of that 
the parameterized wave functions for the deuteron and helium are used. These param- 
eterizations take the spin structure of these nuclei into account. In order to describe 
interactions of the nucleons in the intermediate state, the parameterized A^ A" t-matrix is 
applied that allows us to avoid the problem of convergence which appears at the partial 
wave decomposition at these energies. 

The presented model has been applied to describe differential cross sections at deuteron 
energies of 300 MeV, 493 MeV, and 520 MeV. A reasonable agreement between the data 
and theoretical results has been obtained for the energy equal to 300 MeV. It is shown 
that the contribution of the single scattering term is small at the forward scattering angles 
while inclusion of the rescattering diagram significantly improves the description of the 
experimental data at the scattering angles larger than 30°. The energy dependence of the 
T 2 q has been also obtained at the energy range between 200 MeV and 520 MeV at the 
zero scattering angle. Some improvement of the data description has been received when 
the single scattering term is taken into account. All these results allow us to regard this 
approach as the next step in addition to the one-nucleon-exchange mechanism to solve 
the four-nucleon problem. 
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